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ABSTRACT 

In order to investigate possibilities to measure non-Gaussian signatures of the non-linear iSW 
effect, we study in this work the family of mixed bispectra (r'y-'"*) and trispectra (r'^y"^"^) 
between the integrated Sachs-Wolfe (iSW) temperature perturbation t and the galaxy over- 
density y. We use standard Eulerian perturbation theory restricted to tree level expansion for 
predicting the cosmic matter field. As expected, the spectra are found to decrease in amplitude 
with increasing q. The transition scale between linear domination and the scales, on which 
non-linearities take over, moves to larger scales with increasing number of included iSW 
source fields q. We derive the cumulative signal-to-noise ratios for a combination of Planck 
CMB data and the galaxy sample of a Euclid-liks survey. Including scales down to ^n,ax - 
1000 we find sobering values of cr 0.83 for the mixed bispectrum and cr ^ 0.19 in case of 
the trispectrum for q = 1. For higher values of q the polyspectra (r^y) and (r^y) are found to 
be far below the detection limit. 

Key words: cosmology: large-scale structure, integrated Sachs- Wolfe efi'ect, methods: ana- 
lytical 



1 INTRODUCTION 

The integrated Sachs-Wolfe (iSW) effect is one of the secondary 
anisotropies of the cosmic microwave background (CMB). Time- 
evolving gravitational potentials in the large-scale structure gener- 
ate temperature fluctuations in the CMB ( |Sachs & Wolfe|[l967| l. 
The linear part of this effect is a valuable tool for investigating dark 
energy and non-standard cosmologies since it is sensitive to fluids 
with non-zero equation of state ( Crittenden & Turok|1996 ). For this 
reason its detection is of particular relevance for cosmology and the 
nature of gravity jLue et al.|2004[|Zhang|2006) even though its sig- 
nal strength is very low. 

The linear iSW effect has been measured in such cross- 
correl ation studies ([Bough n et al.||1998| IBoughn & Crittenden 



2004l|Vielva et aLllOOef lMcEwen et al. 2007; Giannant onio et al. 



2008) 1. There are, however, doubts on detection claims formulated 
by |Hemandez-Monteagudol j2010[ > and [Lopez-Corredoira et al.| 
(•2010^ who point out that the iSW-signal seems to be too weak 
on low multipoles below C ~ 10, and that field-to-field fluctuations 
and sampling errors can be important. These facts may correct the 
detection significance to-date to a number of less than two. 

While the linear iSW signal is a large scale effect and becomes 
negligible at angular wave numbers above £ ~ 100, non-linear evo- 
lution of the gravitational potential and leaves signatures on much 
smaller scales, also called Rees-Sciama effect ( |Rees & Sciama| 
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|1968| l and surpasses the linear iSW-effect on these scales. The pos- 
sible signatures of this effect in angular cross spectrum have been 
thoroughly studied analytically (Martinez-Gonzalez et al. 1994) 



|Sanz et al.p996)|Seljak|1996[rSchafer & Bartelmann,2006t . The 



effect increases the total iSW signal by roughly two orders of mag- 
nitude at angular scales around ( ~ 1000 ICooray 2002), before 
gravitational lensing and kinetic Sunyaev-Zel'dovich effect become 
dominant at even smaller scales. However, comparisons of theoret- 
ical studies with numerical simulations showed the Rees-Sciama to 
be negligible in comparison with primary anisotropies on angular 
scales larger than 6 > V (Tuluie & Laguna| ri995 : Seljak ||l996} . 
Also from cross-correlations of the CMB with weak lensing sur- 
veys only a detection significance of ~ 1.5cr from Planck+LSST is 
expected ( [Nishizawa et al.|2008[ l. 

One option to obtain direct signatures of non-Gaussianities 
is the investigation of higher order connected correlators ( Schafer) 
2008 ). In this work we aim to formulate a perturbative approach of 
the mixed iSW-galaxy polyspectra, concentrating on the tree-level 
bispectra and trispectra in flat sky approximation. The unequal rate 
of linear and non-linear evolution at different scales will lead to in- 
teresting sign changes in the spectra, which will also be apparent in 
the non-trivial time evolution of the different source field contribu- 
tions. In addition, we will study the signal-to-noise spectra for mea- 
surements expected from Planck CMB data in cross-correlation 
with observations from a Euclid-like survey assuming unbiased 
measurements with Gaussian noise contributions. We revisit a pre- 
vious estimate of the observability of the iSW-bispectrum jSchafer| 
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|2008 l) correcting an error in the expression for the spectrum of the 
gravitational potential and because of the significantly improved 
signal-to-noise computation, which uses an adaptive Monte-Carlo 
integration scheme l |Hahn|[2005f instead of a binned summation 
over the multipoles. 

The article has the following structure: In Section |2] we will 
lay out the theoretical framework for linear and non-linear structure 
formation (Section [2. 3 1 and Section [2!4| , as well as for the theory 
of higher order correlators of the density field in Section [23] Fur- 
thermore, the main fields of interest, the galaxy number distribu- 
tion (Section [2. 6^ and the iSW-effect (Section [2. 6| l are introduced. 
The mixed bispectra and trispectra are discussed in Section[3] with 
specific studies of their weighting functions (Section [3.4[ ( and their 
time evolution (Sectio n 1 3. 5^ . In Section |4] we present the relevant 
noise sources (Section |4. 1^ , the resulting covariances (Section [4.2^ 
of the polyspectra and finally derive their signal-to-noise ratios 
(Section \4.3\ . Our results are summarized and discussed in Sec- 
tion|5] 

The reference cosmological model used is a spatially flat 
ACDM cosmology with Gaussian adiabatic initial perturbations in 
the cold dark matter density field. The specific parameter choices 
are = 0.25, = 1, (Tg = 0.8, Qb = 0.04 and Hq = 
100 h km/s/Mpc, with h = 0.72. 



2 FOUNDATIONS 

2.1 Dark energy cosmologies 

In spatially flat dark energy cosmologies with the matter density 
parameter Q.^, the Hubble function H(a) = dlna/dt is given by 

HHa) 



-3(l+u.) 



(1) 



with a constant dark energy equation of state parameter w. The 
value w = -1 corresponds to the cosmological constant A. The 
relation between comoving distance x and scale factor a is given 
by 



r' da ^ 
J„ a-H(a) ' 



(2) 



in units of the Hubble distance = c/Hq 



2.2 CDM power spectrum 

The CDM density power spectrum P(k) describes the fluctua- 
tion amplitude of the Gaussian homogeneous density field 5(k), 
((5(k)5*(k')> = (iTrfSoik - k')/'(fc), and the power spectrum of lin- 
ear evolving density fields SlOi) is given by the ansatz 

PL(k)ock'"THk), (3) 

with the transfer function T{k). In low-flni cosmologies T(k) is ap- 
proximated with the fit proposed by |Bardeen et al.|jl986[ >, 

ln(l -1-2.34^) 



r(<?) 



2.34 g 



X [l + 3.89? + (16.1 qf- + (5.46qf + (6.71 qf] (4) 

where the wave number k = qT is rescaled with the shape parame- 
ter r <S ugiyama|I995[ l which assumes corrections due to the baryon 
density fib. 



r = Qmh exp 



-fib 1 + 



(5) 



The spectrum P(k) is normalized to the variance erg on the scale 
R = S Mpc/h, 

(6) 



= ^ J dkk'P(k)wHkR) 



with a Fourier transformed spherical top hat filter function, W(x) = 
3 ji(x)/x. jc(x) is the spherical Bessel function of the first kind of 
order £ jAbramowitz & Stegun|I97^ . 



2.3 Linear structure growth 

The linear homogeneous growth of the relative density perturbation 
field Sl(x, a) is described by the growth function D+{a) 

5L(k,a) = £)+(«) 5L(k, fl = 1), (7) 

which is the solution to the growth equation i Turner & White|I997[ 



[Wang & Steinhardt|I998[|Linder & Jenkins|2003[ >' 

^DAa) +-(3+ ^DAa) = n„(fl)D+(a) . (8) 

da^ a\ dlna / da 2 a- 

The growth equation can be obtained by combining the linearized 
structure formation equations consisting of Poisson equation, Euler 
equation and continuity equation. 

Considering the special case of flat SCDM cosmologies, 
where = 1, the Hubble function scales like H = Hq a"^'^. The 
growing mode solution then gives the very simple growth function 
D+(a) = a, which is even in more complex cosmologies a good 
approximation during the matter domination era. 



2.4 Non-linear structure formation 

In order to describe non-Gaussianities in the density source field 
generated by non-linear evolution a theoretical approximation is 
required. 

We employ non-linear solutions to the density field from 
standard Eulerian perturbation theory ( [Sahni & Coles| |I995[ 
|Bemardeau et al.|2002| >. One expands the density contrast S(k, a) in 
n-th order perturbative contributions (5*"'(k, a), which can be writ- 
ten in terms of the perturbation theory kernels F'"'(ki, ...,k„) and 
the initial linear fields ^*"(k) = (5L(k, a = 1): 



6(k,a) = YjD'lS'-"\k) 

;i=l 

5<">(k) = jd'q,...jd^q„SD(k-qi.„) 

xF<'"(qi,...,q„)5<"(q,)...<5"'(q„) 



(9) 



(10) 



with qi „ = qi + ... + q„. By inserting eqs. (|9jT0j into the 
evolution equations one finds recursion relations for the ker- 
nels F*"'(qi, q„) by combinatorics i Goroff^ et al.||l986 Jain & 
|Bertschinger| 1994) 1. The explicit symmetrized expressions for the 
second order perturbation theory kernels take a very simple and 
intuitive form: 



F<''(ki,k2) 



G<2'(ki,k2) 



1 k, -k, (k 



2 (ki • k2)2 



7 2 kik2 \k2 ki) 7 kjkl 

4 (ki • kj)- 
7 ■ 2 kiki \k2 ' kj ' 7 kjkl 



1 ki •k2 ki 



(11) 



(12) 



where F*'' = G*'' = 1. One can see that mode-coupling to second 
order reaches its maximum when the contributing modes ki and 
k2 are aligned, whereas the kernel vanishes for anti-parallel modes. 
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When in eqn. \\Q) n different modes qi ...q„ contribute to a mode k, 
wave number conservation holds, enforced by the (5d -distribution: 
k = qi + ... + q„. 



2.5 The n-point functions in perturbation theory 

For an analytic expression of the perturbation theory ji-point func- 
tion one has to expand the fields in the correlator. Due to the as- 
sumed Gaussianity of the initial field (5*'* the correlators with an 
even number of fields (5*'' will later simplify to products of initial 
two-point functions Pi^{k), while all uneven contributions vanish: 



(13) 



For simplicity we use in this subsection the notation S„ = <5(k„). 
Simple truncation of the expansion in eqn. l|9j would lead to an in- 
consistent inclusion of powers of the linear power spectrum P\^{k). 
It is more sensitive to take into account all terms up to a certain 
power m in the linear power spectrum, which is equivalent to in- 
cluding terms with initial fields up to powers 2m. 

In this work we exclusively use tree-level perturbation theory, 
i.e. no perturbative terms with wave number integrations are taken 
into account. Following this path, the density bispectrum g^i-''^-''^ 
can be written as 



^k.fek, ^ 2FW(k,,k2)/'L(A:,)/'Lfe)-^cycl.{l, 2, 3). 



(14) 



The non-Gaussian part of the 4-point function is the trispectrum 
j.^ki,k2,k3,k4 convenient to split its tree-level expression up into 
two parts. The first contribution originates from second order per- 
turbation theory. In this case, two of the fields in the correlator have 
been expanded to second order. The expressions in terms of the ini- 
tial power spectra and the second order kernels are of the type 

?<''((ki,k2),(k3,k4)) = ADlPi^{h)Pi^{k^)x (15) 

( FP)(ki3, -ks) F<''(k24, -k4) PL(fcl3) 
+F»>(ki4, -k4) F<2)(k23, -k3) Pl(^14)) • 

The second contribution is due to third order perturbation theory. 
Here, one field is expanded to third order while the other three re- 
main at linear order. For this reason only one perturbation kernel 
appears in the expression for this type of contributions 

fP'(ki,k,,k3,k4) = 6Dl F^'\^,M,\^i)PL(h)Pdk2)PL(h) ■ 

(16) 

With these two functions the connected perturbation theory four- 
point function up to third order in the linear power spectrum P]^{k) 
can be expressed by the following two tree-level contributions 



,ki,k2,k3,k4 



f*''((ki,k2),(k3,k4))-i-all pairs e {1, 2, 3, 4) 



+ ?'''(k|,k2,k3,k4)+ cycl.{l, 2, 3, 4). 



(17) 



The second order and the third order contributions of T^, however, 
have the same time dependence D\. We will see, that this is not 
longer the case for mixed trispectra. 



2.6 Galaxy distribution 

Galaxies form when strong peaks in the density field decouple from 
the Hubble expansion due to self-gravity. These so called proto- 
halos undergo an elliptical collapse ( |Mo et aLl|1997[ |Sheth et al.| 
[200T| |. 

In contrary to the pressure-less dark matter component the 



baryons inside a dark matter halo can loose energy via radiative 
cooling and form stars. Because of the more complex behavior of 
baryons, the relation between the fractional perturbation An/{n) in 
the mean number density of galaxies (n) and the dark matter over- 
density S = lS.pl p is not yet understood. In a very simple way, how- 
ever, the linear relation between the two entities. 



An 



Ap 



(18) 



is a good approximation in most cases and was proposed by 
[Bardeen et al.| l|1986). The bias parameter b can generally depend 
on a number of variables but for simplicity we set the galaxy bias to 
unity throughout this paper, b = I. An established parametrization 
of the redshift distribution n(z) dz of galaxies is 



n(z) dz = no 



exp 



dz with 



1 

"0 



(I) 



(19) 



which was introduced by |Smail et al.| ( [T995^ and will also be used 
in this work. The parameter zo is related to the median redshift of 
the galaxy sample Znicd = 1.406 zo if j8 = 3/2. For Euclid the me- 
dian redshift is Zmed = 0.9. Finally, the F-function ( [Abramowitz &| 
|Stegun|1972^ determines the normalization parameter itq. 



2.7 ISW-efTect 

Due to its expansion our universe had cooled down sufficiently to 
allow the formation of hydrogen atoms at a redshift of z ^ 1089 
( [Spergel et al.|[2003] l. Fluctuations in the gravitational potential 
imposed a shift in the decoupled photons which were emitted in 
the (re)combination process (Sachs- Wolfe effect). This primary 
anisotropy can be observed in the cosmic microwave background 
(CMB) in form of temperature fluctuations AT/Tqmb - 10^ on 
large scales around its mean temperature Tcmb = 2.726 K l |Fixsen| 
[2009 ). 

Besides this, photons are subjected to several other effects on 
their way to us, which lead to secondary anisotropics (reviewed 
by [Aghanim et al.|2008[ l, of which only the most important ones 
are mentioned here: Gravitational lensing ( |Hu||200^ , Compton- 
coUisions with free cluster electrons (Sunyaev-Zeldovich effect, 
[Zeldovich & Sunyaev|[T980) and with electrons in uncoUapsed 
structures (Ostriker-Vishniac effect, Ostriker & Vishniacll986 l and 



gravitational coupling to linear jSachs &~ Wolfe 1967 integrated 
Sachs-Wolfe effect,) and non-linear time-evolving potential wells 
( jRees & Sciama|1968| . Heuristically, the latter two effects originate 
from an unbalance between the photon's blue-shift when entering 
a time varying potential well and the red-shift experienced at the 
exit. 

Assuming a completely transparent space, i.e. vanishing opti- 
cal depth due to Compton scattering, the temperature fluctuations 
t(6) generated by the iSW-effect can be expressed by the line of 
sight integral jSachs & Wolfe|1967[ l 



r(e): 



AT,, 



2 

? Jo 



oa 



(20) 



reaching out to the limit of Newtonian gravity. Using the Poisson 
equation we can write this integral in terms of the dimensionless 
potential = A"' 5lx\ of the density field 6(0x,x)- The n-th per- 
turbative order of the iSW temperature fluctuation r : 
can now be written as 



r(l)4_^(2). 



T<">(6») : 



3n 



, A-' / d Df\ 
dxa'H(a) — \-^\6{0x,x)- (21) 
c Jo Xh V'J" ^ 
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3 MIXED BISPECTRA AND TRISPECTRA 




0.3 

scale factor a 



Figure 1. Time evolution functions Qy'\a) (solid lines) and Q^"\a) (dashed 
lines) as a function of the scale factor a for different perturbative orders 
n = 1,2,3. 



The linear effect {n = \) vanishes identically in matter dominated 
universes = 1, since then D+la is a constant. Therefore, a non- 
zero iSW-signal will be an indicator of a cosmological fluid with 
w Q. After the radiation dominated era it will thus be a valu- 
able tool for investigating dark energy cosmologies. The non-linear 
contributions {n > 2) are now sourced by time derivatives of the 
higher perturbative orders of the gravitational potential. Therefore, 
the Rees-Sciama effect is also present in SCDM-cosmology. 

In order to identify the sources of the effect it is sensible to 
investigate the cross correlation of the iSW amplitude with the line 
of sight projected relative galaxy over-density y = y*'' -I- y*^' -I- ... 



(22) 



JO Ax 

Since we are interested in rather small scales, where non- 
linear effects appear, one can approximate the sphere locally as be- 
ing plane and perform a Fourier transform 



0) 



(23) 



The observable t can be transformed in analogous way. For later 
notational convenience we define the two weighting functions 



bn{z) 



dz 

Ax 



3 fira 



H 



and the time evolution functions to n-th order 



da I a 



(24) 



(25) 



In Fig.[T|the time evolution functions 2y'(a) and (^"^{a) are de- 
picted in dependence on the scale factor a. It is particularly inter- 
esting to observe the different signs. While the galaxy spectra are 
always positive, the iSW contributions change their signs with per- 
turbative order. As we will later observe, also the signs of «-point 
functions will change consequently in the transition from large 
scales, where the linear theory is valid, to small scales, where non- 
linearities start to dominate. For the cross-bispectrum, this effect 
has already been studied ^Nishizawa et al.|2008^ . 



3.1 The density polyspectra 

Regardless of the existence of initial non-Gaussianities in the 
density field ^(k), non-linear structure formation leads to non- 
vanishing three-point and higher order correlators due to quadratic 
terms in the continuity and Euler equation. Since a Gaussian field 
can uniquely be represented by its two-point correlator <f(r) = 
((5(x)(5(x + r)>, multi-point correlators represent a convenient mea- 
sure of evolving non-Gaussianities. The Fourier transforms of these 
2-point and 3-point correlators are the bispectrum g^'-''^*' and the 

, - , ki ,k2.ki,k4 

tnspectrum 

<,5(k,)5(k2)<5(k3)) 

<5(ki)(5(k2)5(k3)^(k4)). 



(2;r)^5D(ki...3)S^ 
(2;r)^'(5D(ki...4) T, 



ki,k2,k3 



(26) 



where the Dirac ^o-function is a result of homogeneity. 



3.2 Limber Projection 

In the flat sky approximation one can use a simplified Limber pro- 
jection jLimber||l953^ to relate the 3-dimensional source spec- 
tra RkbRTHrr^J^ jnrri2.k3.k4 ar,rt„lQr o>^o^tro oh-h-h 



tra b; 



to the angular spectra By 



and 



•XH J 

Ax-,W'y{x)D%(a)T^-''^-'''-'^. 
X 



(27) 



Then, equivalent formulae as in eqn. (|26J apply to these angular 
polyspectra, which are then related to the projected density field 
y{t) with two-dimensional angular wave vectors f,: 



(7(^)7(^2)7(^3 )r(^4)> 



7 



(28) 



where the fields on the sphere with angular directions n, are simply 
decomposed into Fourier harmonics instead of spherical harmonics 



(y(ni)...y(n„)) 



J (27r)2-J (iTiy 
■<y(f,)...7(^«))e'''-"' 



(29) 



Since the region on a sphere around a certain point can for small 
angles be approximated by the tangential plane, this is a good ap- 
proximation for high ^-values. It can generally be transformed to 
the full sky representation with Wigner 3j-symbols ( |Hu|200iy 



3.3 Mixed iSW-galaxy polyspectra 

An equivalent procedure of definitions as in the previous subsection 
can be applied to the iSW-fields T(t). However, due to the uncorre- 
lated noise sources in the iSW and galaxy fields mixed spectra are 
of predominant interest to us. If there exists a chance to securely 
measure the iSW signal it will only work via its cross-correlation 
to the projected galaxy density field y(t) in the cross power spec- 
trum and higher order correlators. 

To allow a compact definition of the mixed spectra we intro- 
duce a doublet field (pj(f) 



7(e) 

T{{) 



(30) 
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Figure 2. Line-of-sight weighting functions W^q^ix) for mixed iSW-galaxy 
trispectra as a function of comoving distance. 



(31) 



Mixed spectra can now be defined in a compact way 
3.4 Weighting functions 

For a mixed n-point function, the product of the n different weight- 
ing functions, is uniquely given by the sum of the field indices q. In 
case of the bispectrum we would define q = ii + 12 + ^3, whereas in 
case of the trispectrum q = ii + i2 + + in- We can therefore define 
a ^-dependent combined weighting function M'^"'(^) 



w?Cr)iv;-''Or) 



(32) 



where w = 3 and n = 4 correspond to the bispectra and trispectra, 
respectively. The different weightings in case of the trispecfi'a are 
depicted in Fig.|2]for different field mixtures q. Despite the weight- 
ings show strong differences in amplitude and sign, common to all 
weightings is a broad peak between 1 and 4 Gpc /i"' due to the 
maximum in the galaxy redshift distribution p{z)- 



3.5 Time evolution 

The time evolution of each linear galaxy field 7(0, k) is given by the 
growth function D+{a}. The n-th non-linear perturbative contribu- 
tions evolve simply with the n-th order of the growth function £)" . 
This is not the case for the iSW field contributions T(a, k). While 
the linear term evolves proportional to d(D+/a)/da the higher or- 
ders can not just be written as the n-th power of the linear growth 
but are proportional to d(D'l/a)/da. 

Due to this fact, different perturbative contributions to mixed 
bispectra and trispectra will in general not have the same time evo- 
lution. In order to obtain a compact notation we introduce the time 
evolution doublet to n-th order Q*"*(a) 



Q*"'(a) 



(33) 



With these time evolution functions Q'"'(a) we are now able to 
write down the general mixed time evolving source fields. For the 



tree-level bispectra we define 



B 



ki,k2,k3 _ 



k3.ki 

■^6 



(34) 



The terms (xh^iT^' are the Poisson factors from the iSW effect. 
In case of the tree-level trispectrum the source will consist of two 
contributions - one originating from second order and third order 
perturbation theory respectively. The time dependent source for the 
trispectra then reads 



ki,k2,k3,k3 



-I- all pairs 6 (1,2,3,4) 



+ Qr,ia)Q\'>{a)Q\'^{a)Q\'>{a) 
-I- cyclic {1, 2, 3, 4) ) . 



,(3)(ki,k2,k3,k4) 



(35) 



Now, the flat sky Limber equations for the mixed angular bispectra 
and trispectra read ( |Hu|200l] l 



'1'2'3 



^''1.^2/3/4 
il'2'3'4 



X 
X 



•XH 1 

dx-W'fHx)Bl 
X 

nXH 1 

dX-Wf>(x)T! 
X 



ki,k2,k3 

'2 '"3 



ki ,k2,k3,k4 
12 13/4 ' 



(36) 



where the source field spectra are evaluated at the 3-dimensional 
wave vectors k, = (/,,i, /,,2, 0). Since the weighting functions are 
slowly varying in comparison to the source field, fluctuations in the 
line-of-sight direction are sineared out by the integrations. There- 
fore, the fields can be assumed as non-fluctuating in this direction 
in the first place. 

While pure spectra are invariant under exchange of wave vec- 
tors {j. 



^aaa 



B. 



B. 



.f4.f1.f2.f3 



,f».f*.f. 



(37) 
How- 



this does not generally hold true for mixed spectra i?^^. 
ever, in general, all spectra are invariant under a simultaneous ex- 
change of wave numbers and field indices 

f1.f2.f3 ^ Df2.f3.f1 ^ Df3.f1.f2 

'2'3'1 '3'1'2 
_ J. f2.f3.f4.f 1 _ 'T.f3.f4.fl.f2 _ J- f4.f 1 .f2.f3 
~ 12/31411 ~ I3'4'1'2 ~ '4'1'2'3 



B 

111213 

y fl .f2.f3.f4 
1112/314 



(38) 



These symmetries are simply caused by the commutation invari- 
ance in the products of source fields in eqn. \26\ . 

3.6 Equilateral bispectra and square trispectra 

To require homogeneity the wave vector arguments have to form 
a triangle, (i + C2 + ti =0, for the bispectrum and a quadrangle, 
(i+(2 + fi + {4 = 0, incase of the trispecfi'um. Thus, due to isotropy 
the scale dependence of the bispectrum is uniquely defined by the 
absolute values of the angular wave vectors f,- 



B 



abc 



B 



la.lb.tc 
abc 



(39) 



The scale dependence of the trispectra, however, can be described 
by the four absolute values of the angular wave vectors ti and one 
diagonal L 



j,ea.tb.ec.ij 

abc 



J. [„.(l,.[c.[j.L 
abc 



(40) 
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Figure 3. Time evolution functions for equilateral mixed iSW-galaxy bis- 
pectra as a function of angular scale C. The value q = ii + ('2 + (3 

defines the mixture of the source fields. While the growth functions of the 
galaxy distribution stay positive to all perturbative orders, the derivatives in 
the iSW evolution functions also introduce negative terms into the evolu- 
tion. 



Figure 4. Time evolution functions for square mixed iSW-galaxy trispectra 
T,.' ■ ' : , as a function of angular scale C. The value a = ii + ij + 13 + u 
defines the mixture of the source fields. The solid Hues depict the second or- 
der perturbative time evolutions S^^uaie- the third order terms Sg'stiuarc ^''^ 
show as dashed lines. While the growth functions of the galaxy distribution 
stay positive to all perturbative orders, the derivatives in the iSW evolution 
functions also introduce negative terms into the evolution. 



This leads to the fact that the source fields of equilateral bispectra 
are symmetric with respect to their field indices and have a uniform 
time evolution Q^^equiia) 



Qq.ei 



\da a j \da a ) 



(9 = 0) 

(9 = 1) 

(9 = 2) 
(9 = 3) 



(41) 

These time evolutions are depicted in Fig. |3] While the growth 
functions of the galaxy distribution stay positive to all perturbative 
orders, the derivatives in the iSW evolution functions also introduce 
negative terms into the evolution. This will later lead to a change 
from correlation to anti-correlation along the line-of-sight. 

Slightly more complex is the time evolution for source fields 
of the square trispectra. Here, the contributions from second order 
perturbation theory evolve still differently compared to the third 
order terms. The second order terms (2Sqi"are('3) read 



1(2) 



6^+(da~a~) 3 ^+ ( do IT ) ( do V ) 

2 +yAa a I \da a / 2 + \(ki a j\<la a ) 
I d« a j \Aa a } 



0) 

1) 



(9 = 2) 
(9 = 3) 
(9 = 4) 



Q(3) 



(q = 


0) 


(9 = 


1) 


(9 = 


2) 


(9 = 


3) 


(9 = 


4) 



and for the third order we obtain the following time evolution 

2^+{da~) 2^+ (da ~o~) (da it) 
4 + yda a j 4 + yila a j\da a J 

(a£L\(as±]^ 

Ida a J\da a ) 

In time evolutions for the different perturbative orders of the mixed 
trispectra are depicted in Fig.|4] One can see that for g = 1,2 they 
evolve identically but differ stronger with increasing number q of 
included iSW fields. 

The equilateral bispectra and square trispectra are depicted in 
Fig. [5] and Fig. |6] As in the power spectra one can observe also 
here the weakness of the iSW signal in comparison to the projected 
galaxy distribution field. This is clearly shown in the decrease of 
the polyspectra with increasing number q of included iSW source 
fields. Once more, the iSW effect shows its nature of being a large 
scale efi"ect. With higher q the slope of the spectra increases in the 
large t region. The physical reason for this is the coupling of the 
iSW effect to the gravitational potential in contrast to the galaxy 
distribution, which couples directly to the density contrast. Math- 
ematically, this fact manifests itself in the appearance of the \lk^ 
factors for the iSW contributions, originating from the inversion of 
the Poisson equation. 

While the linear iSW effect is strictly anti-correlated with re- 
spect to the galaxy density, cross-spectra with this linear signal 
would never show a change in sign in dependence on t. This does 
not hold true any longer for the non-linear iSW efi'ect. The higher 
order contributions now have the opposite sign in their time evolu- 
tion. It is therefore possible that for large I the linear effect domi- 
nates while for small E the non-linear effect determines the sign of 
the correlation. These changes in sign in dependence of t can now 
be observed in Fig. [5] and Fig.|6] for instance at f « 80 in the bis- 
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Figure 5. The absolute values of the mixed equilateral iSW-galaxy bispectra 
^ififi^ as a function of angular scale C are depicted in this figure. The value 
1 = h + h + '3 defines the mixture of the source fields. 




Figure 6. The absolute values of the mixed square iSW-galaxy trispectra 

T. as a function of angular scale ( are shown in this plot. The value 
II 12/314 " 

q = il + (2 + !3 + h defines the mixture of the source fields. 



pectrum (ry-). This behavior has also been observed and studied 
in the CMB-weak-lensing cross spectrum ( [Nishizavya et al.|2008^ . 



4 DETECTABILITY 
4.1 Sources of noise 

The step from a good theoretical framework to an analysis of real 
data or to an estimation of the realistically accessible inforination 
content encompasses the description of all relevant effects influ- 
encing the measured data. Only then, one will be able to make 
statements about a physical process and the likelihood of its actual 
measurement. For the evaluation of the covariances of the bispectra 
and trispectra, the two-point function will be needed. In the same 
notation as the higher order spectra their are defined as 

<</)„ (e,)ip„({2)) = (2;r)2<5D(/i + {2) (42) 

In the case at hand the actual theoretically expected iSW signal in 
our fiducial cosmological model is superposed to the primary CMB 



N 


An 


/sky ZO 


h n 


3.0 X 10' 


In 


0.5 0.64 


1.0 4.7x10'* 



Table 1. The Properties of the Euclid galaxy survey are listed in this table: 
total number N of objects, solid angle Af2 covered (in radians), sky fraction 
/sicy, redshift parameter zo, galaxy bias b and density per unit steradian n. 



fluctuations. Its relative amplitude reaches from 10% on very large 
scales to an negligible fraction of the signal for scales smaller than 
E x 200. Furthermore, the detected CMB signal is subjected to 
instrumental noise cr^ and a Gaussian beam p{£). 

Assuming that the noise sources of the galaxy counts are mu- 
tually uncorrelated, the pure galaxy-galaxy spectra are solely sub- 
jected to a Poissonian noise term . 

The cross-spectra between the two fields will be free of noise, 
since the noise sources of the single fields are uncorrelated. 

Now we can relate the measured spectra C^'^ to the theoretical 
spectra : 



t 

00 

t 
11 



c 



C„n + n 



'-Ol 



C + C 



+ oii3-\n. 



(43) 



The contributions in detail are: 

(i) As the Fourier transform of the Gaussian beam one ob- 
tains p-^it) = exp(Ae^{(£ + 1)). We use AO = 7.1 arcmin, which 
corresponds to the v = 143 GHz channels closest to the CMB 
emission maximum. For the conversion of w^' = T'cmb""? 'he 
noise amplitude in the dimensionless temperature perturbation t 
with wt = (0.01 fiKf ( |Zaldarriaga et al.|l997| we use the value 
TcMB = 2.725 K for the CMB temperature. 

(ii) Furthermore, a CMB temperature power spectrum 
CcmbC^) was generated, which was equally scaled with the CMB 
temperature Tcmb = 2.725 K, with the Code for Anisotropics in 
the Microwave Background (CAMB, [Lewis et aL]|2000| l for the 
fiducial ACDM cosmology. The noise contribution from the CMB- 
spectrum C[.^g represents the main challenge in the observation of 
the iSW bispectra and trispectra. It provides high values for the 
covariance at low multipoles C, and it by far dominates Coo(i), 
CcMB s> Coo(/) on the angular scales considered. The orders of 
magnitude for the different contributions of the linear estimator Cj , 
are depicted in Fig|7] One can see, that even at large angular scales 
C the pure iSW signal C{j is still more than one order of magnitude 
weaker than the signal from primordial fluctuations C^j^^- 

(iii) The inverse number density n of objects per unit stera- 
dian determines the Poissonian noise term in the galaxy count. In 
Table |4.l] the properties of the main galaxy sample as it would be 
expected from Euclid are summarized. Major advantages lie in the 
large sky coverage and the high number of observed objects. Here, 
we assumed a non-evolving galaxy bias for simplicity. 



4.2 Covariances 

In the case of Gaussian noise the observed and estimated bispectra 



fi.'l''^''^ and the trispectra f 

'1'2'3 '1'2'3'4 



are unbiased estimates of the 



true bispectra B I' - 

'1'2'3 



and trispectra y 'i-'^i's-'"'! 

^ ll'2'3'4 



iHu 



2001 1, 



'1'2'3'4 



h'2'y 

■T^^l'^2'^3'''4 
'1'2'3'4 



(44) 
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Figure 7. Constituents of the measured angular CMB spectrum Cj j . De- 
picted are the total signal Cf j (solid line), the contribution from primordial 
fluctuations C^q^^, the iSW-elfect Cj j and the instrumental noise cr^ p^^(t), 
which is fortunately sub-dominant at the large scales of interest. 



q 





1 


2 


jii) 


87.8 


0.828 


4.43 ■ 10-3 


^1 


21.7 


0.19 


1.42- 10-3 



Table 2. Cumulative signal-to-noise ratios Z^"' for measurements of the bis- 
pectra {t^y'^i) and the trispectra (i^y^'i), g = 0, 1,2, for Planck CMB 
data in cross-correlation with Euclid-liks survey, up to a resolution limit 
^max — 1 03 starting from a minimum angular wave number of ^min — 10. 



trispectra in tlie subspace £i , E\ < ^2, i'l < ^3. ^3 < A> ^4 amounts to 



Cov-' If f'- 'i-'-''' 

' "hrri ' n' h' r' rl' 



^nn> ^hh' ^^J' 



det C 'i det C ^2 det C det C 
X Soiti - t\ ) 5d(^2 - t'l) <5d(^3 - ^3) '^D«'4 - O ■ (49) 
For an observation covering the sky with a fraction of /s^y the co- 
variances scale like /.^J,. The anti-correlation in the cross-spectra 
Cqj will not change the sign of the covariances, since in each of the 
products an even number of these mixed spectra appears. 



This is in contrast to the spectra C\^~^, which were discussed in the 
previous subsection. The covariances of the estimators of the bis- 
pectra and trispectra are defined as 



Cov S 



f1.f2.f3 h'W'z 



11(2/3 '1'2'3 / \^ M'2'3 '1'2'3 // ' 



(112(3,14 

\\ '1'2'3.'4 ill2'3.'4 ) \ 'ihi-iM '1'2"3.'4 // ' ^ ^ 

In a Gaussian approximation, which we are using here, any covari- 
ances can be expressed as a sum of products of two-point func- 
tions using Wick's theorem. While for pure covariances only the 
respective power spectra appear in this expansion, in our case of 
mixed covariances the products are formed from the estimators of 
the cross-correlation Cgj and the two auto-correlations C^q and C[j . 

In case of the bispectra with mutually unequal angular wave 
numbers # for j + k the covariance can be written as a sum 
over terms which are cubic in the spectra C\^._^ 



Cov B 



f 1/2/3 

abc 



~ f f f 
a'b'r 



+ p&Tm{l\J'2J\). (46) 

On the subspace £i,t\ < 1^,1'^ < t^t'^ only the first term is non- 
vanishing. This block-diagonal matrix can now be inverted to 



aa bb cc 



[ ' \ detC^idetC'2detC'3 

X feC^i - l\ ) 6u{t2 - {'2) 6u(h -t'i), (47) 
with the adjoint matrix C'J^,, 



^11 ^01 



C"=\ _ V . (48) 

\ "^01 '-00 / 

In analogy to the bispectrum case, the inverse covariance of the 



4.3 Signal-to-noise ratios 

The signal-to-noise ratio E*3) for the simultaneous measurements of 
the all pure and mixed bispectra {r^y^^'') and S^*' for the all mixed 
and pure trispectra {r^y*^''), where all field indices are summed 
over, would imply a thorough derivation of all cross-correlations 
between different field mixtures. Here, we are rather interested in 
the individual signal-to-noise ratios of certain field configurations. 
If one reduces the data to a mixed configuration q, only the cases 
^ = and q = \ provide measurement in and above the detection 
limit. For ^ = we obtain 



(^(3,V 







; 



,/sky 
4^ 

/sky 

9,n^ 



! 



[^000 ) 



d^^i d-Z'j d^fs 



y 0000 / 
^^'-oo'-oo'-'oo'-oo 



(50) 



However, since we are aiming for iSW detections, the more inter- 
esting case is g = 1. The signal-to-noise ratio then splits up into 
two contributions, 

/sky 
4^ 



K-^ooi ) "^1 1*^1 1 '-00 



''1/2/3/^^1 ^2/^ ^3 D^l'^2/3 



V) 



/sky 
8^ 



J" drtiA^Cid^ti (detC'' 

1 /'y,fl/2.f3/4\^ /=lfl/Sf2/=;f3^ 

£ yoool I '-'u '^ii 



'"ll *^oi "-oi -"oio 



.detC' 



00 



Tifl/2/3/4/^ fl f'hf'hf'Uj- f 1/2/3/4 
'0001 '-'U "-11 ""Ol "^01 ^ 0010 



(51) 



A detailed calculation of the signal-to noise expressions can be 
found in Section |A] The inverse covariances of the polyspectra 
will always remain positive, since always an even number of anti- 
correlating cross-spectra will appear in its expression. However, the 
mixed field contributions can in general become negative. 

The cumulative signal-to-noise ratios E^"' for the mixed bis- 
pectra Bq and the mixed trispectra are depicted in Fig.[8]for the 
pure galaxy spectra and spectra with up to two iSW source fields 
included, ^ = 0, 1, 2. 
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■7 I I 
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Figure 8. Cumulative signal-to-noise ratios Z^"' for measurements of tiie 
bispectra (riy^^i) (solid lines) and the trispectra {r^y^^i) (dashed lines), 
I? = 0, 1, for Planck CMB data in cross-correlation with EucUd-\\ks survey, 
up to a resolution limit f„,ax = 10^ starting from a minimum angular wave 
number of — 

10. 



The qualitative behavior of the cumulative signal-to-noise 
curves are again determined by the individual signal strengths of 
the two source fields y and t. The strong fluctuation of the galaxy 
distribution y even on smaU scales leads to a considerable increase 
of Z for large £ and small q. In contrast to this the iSW-efFect is 
a large scale effect and therefore increases the slope in the small 
£ limit of the spectrum. It does not contribute significant signal 
strength above values of C^.^^ > 300, for this reason the signal-to- 
noise curves flatten off in this region of the spectrum for q = 1,2. 
The wider spread between different values of q for the trispectrum 
in contrast to the bispectrum is due to the higher power of source 
fields. 

Quantitatively, higher values of q lead to smaller significance 
in the signal. Included were contributions starting from large an- 
gular scales £^j„ = 10 up to smallest scales measurable in the 
Planck survey ^„,ax = 10'. The pure galaxy poly spectra (y^) and 
(y*) can both be measured with a detection significance of » 3cr, 

= 87.8 and Zj,'*' = 21.7. Including only one iSW source 
field reduces the signal down to the noise level. While the bispec- 
trum (ry^) reaches a signal-to-noise ratio of 0.82, the value for the 
trispectrum (ry^) reaches a maximum of 0.19. Combining mea- 
surements of the q = 1 bi- and trispectra would therefore be able 
to contribute a maximum signal-to-noise contribution of E k 0.84. 
Unfortunately, this is - taken on its own - still a very poor measure- 
ment significance. However, it could be used as an additional signal 
source to the strongest iSW signal from the cross spectrum (ry). 

For the higher values of q only the case of two iSW source 
fields q = 2 is plotted in Fig. [8] Both for the bispectrum as well 
as for the trispectra the signal-to-noise ratios are negligible with 
maximum values of 4.43 • 10"' and 1.42 • 10"' respectively. 

One can obtain a grasp of the differential contributions of the 
signal-to-noise ratios with respect to angular scale d, if one studies 
the quantity 

(^5,(3) Y ^ B[f '' Jcoy-Hb'/-') 

(^s^U)' - r,''^-''^-^'^Cov-'(rr''^| (52) 




10 10 10 

I 

Figure 9. The differential contributions of the equilateral bispectra to the 
signal-to-noise ratios in dependence on angular wave number C are depicted 
here for different source field mixtures ^ = 0, 1, 2 (solid lines). For q = 1,2 
the contributions from cross-correlations are also shown (dotted lines), as 
they appear for ^ = 1 in the second term of the second line in eqn. JsTJ. One 
can observe the increasing amplitude of the baryonic acoustic oscillations 
for larger q. Also the change in sign can be studied due to the transition 
from linear dominated to non-linear dominated scales. 



In Fig. [9] this differential contribution of equilateral bispectra in 
dependence on £ are depicted for different source field mixtures q. 
The differential contributions of the square trispectra behave qual- 
itatively analogous. Also the change in sign can be studied due 
to the transition from linear dominated to non-linear dominated 
scales. For q = 1,2 the contributions from cross-correlations are 
also shown (dotted lines), as they appear for 17 = I in the second 
term of the second line in eqn. \5\\ . As one can see, these terms are 
subdominant and can be neglected in our case. 

One can observe the increasing amplitude of the baryonic 
acoustic oscillations for larger q, which originate from C[.^^ . The 
falling slopes of the BAO features in the covariance lead to small 
plateaus in the differential contributions for larger t. Since in these 
regions the signal decreases more gently than the covariance, one 
obtains a local increase of signal-to-noise. However, this effect can 
hardly be observed in Fig. [8] 



5 SUMMARY 

The objective of this work is a study of the detectability of non- 
Gaussian signatures in non-linear iSW-effect. Besides the mixed 
bispectra of the form {■fy^^''), ^ = 0, 1, 2, between the galaxy dis- 
tribution y and the iSW temperature perturbation r we also calcu- 
late for the first time the mixed trispectra of the analogous form 
{fy*^''). Both types of spectra were consistently derived in tree- 
level perturbation theory in Newtonian gravity. This implies for 
the bispectra perturbative corrections to second order and for the 
mixed trispectra contributions from second and third order terms. 
Furthermore, we investigated the time evolution of these individual 
3-dimensional source terms, which are in general very diverse. For 
this reason, the time evolution and the configuration dependence 
of a specific class of spectra, equilateral bispectra and the square 
trispectra, were studied. Finally, the achievable signal-to-noise ra- 
tios were derived for measurements cross-correlating Planck data 
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and a galaxy sample, as it would be expected from a wide angle 
survey as Euclid. 

(i) The linear iSW-effect has the time dependence d(£)+/a)/da, 
which makes it sensitive to dark energy but vanishes in SCDM- 
models with Sl^ = 1 and /)+(«) = a. In contrast to this, the non- 
linear contributions to the iSW signal are sensitive to derivatives of 
higher powers of D+{a), namely d(D^/a)/da for second order per- 
turbation theory and d(Dl/a)/da for third order contributions. For 
this reason, the effect does not vanish in matter-dominated epochs. 

(ii) The covariances of the measurements were derived in a 
Gaussian approximation. For the CMB observation the intrinsic 
CMB fluctuations and instrumental noise in form of the pixel noise 
and a Gaussian beam were considered as noise sources. A Pois- 
sonian noise term was added to the galaxy distribution signal. For 
simplicity the fluctuations of the dark matter density and galaxy 
number density were related to each other by a constant linear bi- 
asing model. 

(iii) In the mixed bispectra and trispectra the configuration and 
scale dependence represent the different correlation lengths of the 
gravitational potential and the density field. Since the specific per- 
turbative corrections dominate on different scales, the mixed spec- 
tra change their sign at certain values of £. In case of the bispectra 
one can observe the transition from linear domination to non-linear 
domination move to larger and larger scales with increasing num- 
ber of included iSW source fields q. 

(iv) We derived the cumulative signal-to-noise ratios "Zf^ for 
the measurements of mixed bispectra {t^j^^''), and Xjj'** for the 
mixed trispectra of the form (r'y'*"'), with a Gaussian approxi- 
mation to the covariance. The integration were performed numer- 
ically using Monte Carlo integration techniques from the multidi- 
mensional numerical integration library CUBA (Hahn 2005). For 
both spectra the initial CMB fluctuations are the most important 
noise source, which makes it difficult to observe the signals. We as- 
sumed a cross-correlation of Planck data with a Euclid-\ike. galaxy 
sample starting from angular scales of C^^^ = 10 up to a resolution 
of ^raax = 10'. The only spectra reaching the order of magnitude 
of the noise level are the bispectra and trispectra in the config- 
uration (ry""'). We found the numerical signal-to-noise ratios of 
"Zf^ = 0.828 for the bispectrum and ^ = 0.19 for the trispectrum 
and conclude, that non-Gaussian signatures of the iSW-effect are 
too weak to be detected. At the same time, these small signal-to- 
noise ratios suggest that non-Gaussianities in the CMB generated 
by the iSW-effect are small enough so that they do not interfere 
with the estimation of the inflationary non-Gaussianity parameter 
/nl from the bispectrum (t'> and of the two parameters gNL and 
TfjL from the trispectrum (t"*). 
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APPENDIX A: ANALYTICAL DETAILS OF 
SIGNAL-TO-NOISE RATIOS 

The squared signal-to-noise ratio is given by the between a 
detection and its zero hypothesis. In the course of its calculation, 
one has to ensure that no redundant information is taken into ac- 
count. We present the calculation for the bispectra only, since it 
follows the same argumentation in the case of the trispectra. Ne- 
glecting redundancy due to any symmetries all mixed and pure bis- 
pectra would account for &x^' contribution of 



X 



/sky 



B, 



fi.h 
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where also the sum over all field indices is implied. However, the 
integrand is symmetric in any simultaneous pairwise permutation 
of with (t„,,i„,) and likewise of (f'„,i'„) with (^'j,, This 

type of redundancy can be avoided by constraining the integration 
volumes to < £2 < ^3 and < l'^ < £y Furthermore, the sum 
over field indices may lead to more redundancy, which we encode 

at this point into a multiplicity factor i'-j'^-j. Now, the signal-to- 
noise ratio can be written as 
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In this subspace the covariance matrix can be inverted, as it was 
shown in Section [4~2] Substituting eqn. into eqn. |A2j, we find 
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If one uses 
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one can combine the quadratic terms to one, which is integrated 
over the full f 1 2,3-volume. This can be done for the mixed terms in 
analogy and one is left with the following expression for the signal- 
to-noise ratio, 

(Sf )' = ^ J d^/'i A^li (2 detC'' detC'^ detC'^)^' 
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Following the analog path of argumentation one finds the signal-to- 
noise expressions for the trispectra to be 
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Also the expressions for higher values of q can now be deduced 
with the same techniques. 

This paper has been typeset from a Tj3f/ KTe^ file prepared by the 
author. 



where from now on = —t\ — I2 is implied, if is not inte- 
grated over. If one is now interested in the signal-to-noise ratio 
of particular field mixture, i.e. data with a fixed field configura- 
tion q = ii + ii + h, one can further simplify the expression. For 
pure galaxy contributions, ^ = g' = 0, we can neglect the cross- 
correlation , i.e. Coi = 0, and no redundancy due to field summation 
occurs, igJJJ = 1. One obtains the well-known case ( |Hu|2001| l 
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where the symmetry in the integrand was used to obtain an inte- 
gration over full f-space in combination with the factor 1/6. If 
all mixed spectra with one iSW field are taken into account, i.e. 
q = q' = I, one obtains 9 different contributions due to the field 
index summation. Three contributions are quadratic in identical 
bispectra, ((i,i2,!3) = 12,12), and have multiplicity one. The 
remaining mixed contributions have multiplicity 2, since the inte- 
grand in eqn. |A3| is symmetric under exchange of the primed and 
unprimed index sets {i\,i2, 13) and (ip ij, i'^) as a whole. Therefore 
the multiplicities are 
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